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MAXIMUM LIKELIHOOD DRIFT ESTIMATION FOR THE 
MIXING OF TWO FRACTIONAL BROWNIAN MOTIONS 

Y. MISHURA 


Abstract. We construct the maximum likelihood estimator (MLE) of the un¬ 
known drift parameter 0 S M in the linear model 

= + + t€ [0,r], 

where B^^ and B^'^ are two independent fractional Brownian motions with Hurst 
indices ^ < Hi < H 2 < 1. The formula for MLE is based on the solution of the 
integral equation with weak polar kernel. 


1. Introduction. The elements of stochastic calculus for fBm 
Consider the continuous-time linear model 

X{t) = et + aB^^{t) + t e [0,T], 

where B^^ and B^^ are two independent fractional Brownian motions with Hurst 
indices ^ < R2 < 1, t > 0. Given the sample path of X on [0, T] it is required 

to estimate the unknown drift parameter 0 G M. In the case when Hi = ^ the 
problem was solved in where the authors develop the basic tools for analysis of the 
mixed fBm based on the hltering theory of Gaussian processes. They considered the 
linear regression setting and demonstrated how the maximum likelihood estimator 
can be constructed and studied in the large sample asymptotic regime. 

As the preliminaries, we give some basic facts about stochastic calculus for frac¬ 
tional Brownian motion (for more details see |2], j5], [7], j8], [IQ]). Parameter es¬ 
timation for different models with long memory was studied, among others, in [1], 
0. 0 . E]. Let (G, P) be a complete probability space supporting all stochastic 
processes considered in what follows. Introduce > O}, an adapted frac¬ 

tional Brownian motion with Hurst parameter H G (1/2,1), i.e. a centered Gaussian 
process with the covariance function 

E [B«(«)B«(s)] = i («“ + s"" -\t- s|“) . 

For a (possibly inhnite) interval [0,T], denote by L^[0,T] the completion of the 
space of simple functions /: [0,T] —)• M with respect to the scalar product 

if, 9)h ■= [ [ fit)gis) \t - dsdt, 

Jo Jo 

where an = H{2H — 1). (It is worth to mention that this completion contains not 
only classical functions, but also some distributions.) For a step function of the form 

n—1 

k=0 

1 
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where {to < < • • • < tn } C [0,T], the integral I ^{ f ) of / with respect to is 

dehned by 


pT n-l 

= (i?^(4+i) - 5^(4)) . 

fc =0 

It can be verified that maps isomorphically the space of step functions on [0,T] 
with the scalar product (•, ■)into hence, can be extended to L^[0,T]. 

Now dehne a square integrable kernel 

Knit, s) = J {u — 

1 

where ( 3 h = { b{h-i/2 2-2H) ) where B{-, •) is beta function. The map 

(4|f/)(s) = f{t)dtK{t, s)dt = /(t)t^"^/^(t - s)^~^/‘^dt 

is an isometry between the space of step functions and and can be extended to a 
Hilbert space isomorphism between L|^[0,T] and L^[0,T]. This implies that the 
process 

W(t) = I« ((A-J)-‘l|„,,|) 

is a standard Wiener process on [0,T], moreover, for any / G L|^[0,T], 

(1.1) J^(/) = l\Ky){s)dW{s). 

Jo 

In particular, putting in the last formula / = l[o,t], one gets the following well-known 
representation of B^: 

B^{t)= f Knit, s)dW{s). 

Jo 

Finally, we define the so-called fundamental martingale, or Molchan martingale , 
for B^. In this order, introduce the kernel s) = (f — and consider 

square-integrable Gaussian martingale 

M^{t)= flH{t,s)dB^{s)=jH f s^/^-^dW{s), 

Jo Jo 

with 'jH = (2hf(| — i^)r(3/2 — H)^T{H + |)r(3 — the last equality is due 

to fll.ip . 

The paper is organized as follows. In Section 2 we reduce the main problem to 
the solution of the integral equation with the weak polar kernel and establish the 
existence-uniqueness result for this equation. Section 3 contains an auxiliary result 
concerning the existence and uniqueness of the solution of the corresponding integral 
equation of the 1st kind. We prove this fact directly, by constructing the unique 
solution. 
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2. Main problem 


Now, let \ < Hi < H 2 < t > o|, i = 1,2, be two processes 

defined on the space (12, i?, (5^)t) and Pq be a probability measnre under which B^^ 
and -8^2 are independent, B^"^ is a fractional Brownian motion with Hurst parameter 
i/ 2 , and is a fractional Brownian motion with Hurst parameter Hi and with 
drift i.e., 

aB^^{t) = et + (jB^^{t). 

The probability measure Pq corresponds to the case when 6* = 0. Our main 
problem is the construction of maximum likelihood estimator for 0 G M by the 
observations of the process Z{t) = 9t + aB^^{t) + B^^{t), t G [0,T]. However, the 
form of the process Z (two fBm’s with different Hurst indices) does not allow to 
construct the estimator immediately. To simplify the construction, we apply to Z 
the linear transformation of the following form: 


( 2 . 1 ) 


Y{tj = fl„,{t,s)dZ{s) = - Fi, 2 - + aM"'{t) 

+ [ lHi{t,s)dB^^{s). 


this process is preferable since it involves Gaussian martingale . 


Lemma 2.1. The linear transformation fl2.ip is correctly defined. 

Proof. It is sufficient to establish the existence of the integral Ih^ {t, s)dB^^{s) for 
any t G [0, T], But we have that for any u, s G [0, t] 


therefore 



llliit^ s'^lHiif., u') \u S 


\2H2-2 


dsdu 




2Hi-2 


dsdu 
P|M^i(f)|2 


^ ^2H2-AH,+2 

2-2H 


whence the proof follows. 


□ 


As it was mentioned, process Y is more convenient to deal with since it involves 
martingale with a drift. Furthermore, it follows from the next result that processes 
Z and Y are observed simultaneously, so, we can reduce the original problem to the 
equivalent problem of the construction of maximum likelihood estimator of 0 G M 
basing on the linear transformation Y. 

Lemma 2.2. Processes Z and Y are observed simultaneously. 

Proof. Taking into account fl2.ip . it is enough to present Z via Y. But it follows 
from fl2.ip . from Fubini theorem for integrals w.r.t fBm (Theorem 2.6.5 jS]), and 
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from elementary integral transformations, that for any t G [0,T] 


{t — 2 J lHi{s,u)dZ{u)ds = I M2 

13 ' 


h-Hi 


{t — 2(s —■u)2 ^^dsdZ{u) 


= b{h^--,--H^) / u^-^^dZ{u)= / {t-s)^^-W{s)ds 


whence 






-1 


t rt 



{u - s)^^-^l‘^u^^-^/^dudY{s), 


0 Js 


and the proof follows. 


□ 


Denote for simplicity Bhi '■= — Hi, | — Now the main problem can 

be formulated as follows. Let ^ < Hi < H 2 < 1, = M^^{t),X 2 {t) : = 

Jq Iniit, s)dB^^{s), t > 0}, i = 1,2, be two processes defined on the space 

and Pg be a probability measure under which Xi and X 2 are independent, B^^ is a 

fractional Brownian motion with Hurst parameter H 2 , and Xi is a martingale with 

2 

square characteristics {Xi){t) = 2 - 24 / 1 and with drift ^^ 2 - 24 / 1 ^ 

Xiit) = = ^^^2 -2441 ^ 

a 

Also, denote Xi{t) = Our main problem is the construction of maximum 

likelihood estimator for 0 G M by the observations of the process 

Y{t) = + aXiit) + X2{t). 


Note that under measure Pg the process 


W{t) := W{t) + 


<^7h. (i - j 


is a Wiener process with drift. Denote 6hi = 

By Girsanov’s theorem and independence of Xi and X 2 , 


dPg 

dPn 


= exp 95_ 


544i 


52 


-44i 


dW{s) 




= exp 96hiXi{T) - 




4(1 - H 


4{1-Hiy 

_ rjn2 — 2Hi 


~\2—2H\ 




For technical simplicity, we put a = 1 in what follows. Note that, similarly to the 
linear combination of Wiener process and fBm, considered in |1], this derivative is 
not the likelihood for the problem at hand, because it is not measurable with respect 
to the observed a-algebra 


:= cj{Y{t),te [0,T]} = := a{X{t),t G [0,T]}, 


where X{t) = Xi{t) + X 2 {t). 

We shall proceed as in |1]: let jig be the probability measure induced by Y on 
the space of continuous functions with the supremum topology under probability Pg. 
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Then for any measurable set A jUg(A) = ^(x)/j,odx, where <h(x) is such measurable 


functional 


that ‘h(X) = Eq ^Tj- The latter means that /iq for any 6^ G M. 

Taking into account that Xi = Xi under Pq and the fact that the vector process 
{Xi,X) is Gaussian, we get that the corresponding likelihood function is given by 


( 2 . 2 ) 


LT{X,e) = = Eo(exp [e6H,X,{T) - 


02^2 


Hi 


= exp\e6HiEo{Xi{T)\d^) + 




4(1 - H,) 

rj-i2—2Hi 


T 


2-2Hi 


d 


2 V ' ' 2-2Hi. 

where V{t) = Eo{X^{t) - Eo{X^{t)\^f))\ t G [0,T]. 

Thus, we arrive at the following problem; to hnd the projection PxXi{T) of 
Xi{T) onto {X(t) = Xi(t) + X 2 (t),t G [0,T]}. We recall from Section [1] that 

W,{t)= f t = l,2, 

Jo 

are standard Wiener processes, which are obviously independent. Also from Sec¬ 
tion [T] we have 


(2,3) 


Xdt)=~tH, / B«Ht)= / K„,{t,s)dW 2 (s). 


Then, using fll.ip . we can write 


X2{t)= / KHiMt^^)dW2{s), 


where 


(2.4) KHi,H 2 it, s) = /3 h 2S^^‘^ j — ^^'^du. 

Similarly to fll.ip . we have for / G L^J0,T] 


(2.5) 

where 


f{s)dX2(s) = / (A'l,. „J)(s)<iW,(s), 


(A'L,«,/)(«)= / fmKH,Mt,s)dt 


The projection of Wi(T) onto {X{t),t G [0, T]} is a centered X-measurable Gauss¬ 
ian random variable, therefore, it has a form 

rT 


PxXYT)= f hT{t)dX{t) 

Jo 


with hr G L|^j0,T]. Note that still can be a distribution. This projection for 
all u G [0,T] must satisfy 

(2.6) E [X{u)PxXYT)] = E [W(n)Xi(T)]. 

Using fl2.6p together with independency of Xi and X 2 , we arrive at 


(2.7) 


Xi{u) / hT{t)dXYt) + X 2 {u) / hT{t)dX 2 {t) 


= £[XYu)XYT)]=eHiU 


2-2Hi 
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where Eh = i'h/— 2H). 
From fl2.3p - fl2.7p we get 


( 2 . 8 ) EhiU^ ^=Thi ^ t ( s)s / hT{s)rH^^H2{s,u)ds, 


where 


fsAu 


rHi,H 2 {s, u) = 


dsKHi,H 2 (s, 'v)Khi,H 2 {u, v)dv. 


This kernel can be written alternatively as rHj^^H 2 (ty^) = dtRHi,H 2 ity where 

ptAs 


RHi,H2(t, 5 ) = 

rt rs 


— OiH2 



RHi,H2(ty u)Khi,H 2 {s,u)du = E[X 2 {t)X 2 {s)] 

2—-f^i L, \2H2—2 




0 Jo 


Differentiating fl2.8p with respect to u, we arrive to 


(2.9) 
where 

( 2 . 10 ) 


Th,“' ™ = 7H,ftT('“)M' + / hT(s)k(s,u)ds, 


l-2Hi 


psAu 


k{s,u) = durHi,H 2 {s,u) = / dsKH^^H 2 is,v)duKHi,H 2 {u,v)dv. 


Theorem 2.3. Let H 2 — Hi > Then there exists a sequence T„ —)■ cx) such 
that integral equation fl2.9p has unique solution hj-^ on any interval [0,T„] and 

^ -h2[0,n]. 


Proof. We denote Chi,H 2 constants which values are not so important; their values 
can change from line to line. At hrst, we can apply the changing of variables u = 
s + {t — s)z to f|2.4p and transform the kernel KHi,H 2 {tj s) from fl2.4p to the following 
form; 

( 2 . 11 ) 

/ (1 - Z)5-">(S + ((- s)z)''“-»>z"''-}&. 

Jo 

Then we can differentiate fl2.1ip and after inverse changing of variables we get that 


( 2 . 12 ) 


dtKH,,H2itH) = {H2-Hi)(^ 


K 


Hi,H2 


(t,s) 


t 


+(3h2S^ 


h-H2 


t — s 


(t-r) 




r — s 


,H2-h 


^dr 


Further, we have the following bound for kernel Khi,H 2 {'Ii -s) on the interval [0,T] : 

(2.13) 0 < < ^H2b{^ - Hi,H2 - 

and it follows from fl2.12p and fl2.13p that 

(2.14) 

0 < d,K„, nJt,s) < - Hi,H2 - l)i"’""‘(l “ s)""""'"' 

+b(5 - HuH2 + i)««>-«--‘(( - s)"-"‘) < 
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Now we can substitute the bound from fl2.14j) into fl2.inj) and get that 

(2.15) 

psAu 

0 < k{s,u) < 

Jo 

Let, for example, s < u. Note that 

{U - < („ _ _ s)l-2-f^i 

Then it follows from fl2.15p that 
0 < k{s,u) < 

(2.16) 


X 2772 (.y _ 

Jo 


H 2 +H 1-2 


{s-v) 


H 2 -H 1 -I 


dv. 


In order to bound the integral in the right-hand side of fl2.16p . we apply the hrst 
statement from lemma 2.2 |2], according to which for /x, z/ > 0, c > 1 

[ t^-\i - tY-\c - = B{fi, u)c-yc - 1 )-^. 


Therefore, with /i = 2 — 2 H 2 , u = H 2 — Hi and c = - 




(2,17) 




(1 - !>) 


H 2 -H 1 -I 


dv 


1 /-U\I7i-H 2 /-y \ 2772-2 

= B{2 - 2 H 2 , H 2 - Hi)s-^ (-) (- - 1' 

_ (J gl—772 —77i.,,77i—772 7 .,, ^\2H2—2 


-771,772*- " ^Y^-sY 

and it follows from fl2.16p and (12.17p that for s < u 

(2.18) 0 < k{s,u) < Ch,,H 2 S^-^^Yu - 
Evidently, for u < s 

(2.19) 0 < k{s,u) < Ch^,H 2 U^~^^Ys - 
Now we rewrite equation (12.91) in the equivalent form 


( 2 . 20 ) 7771 "^^ = 7 hiY{u)u^ + I h{s)s^ ^k{s,u)ds, 

Jo 


,i-77i 


i-77i^77i-i ,77i-J 


or 


( 2 . 21 ) 


I hT{s)ki{s,u)ds, 


where ki{s,u) = 2k{s,u), hT{u) = h{u)u2 and it follows from (I2.18P 

and (12.191) that for s < u 


ki{s,u) < Chi,H2'^^^ 252 ^^{u — s) 


'\2i^2—— 1 


kl{s,u) < Chi,H2^^^ 2-U2 'f^l(s_.y) 


\2H2—2Hi — 1 


and for M < s 
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Therefore, taking into account that for H 2 —H 1 > | we have that AH2—4:Hi—2 > —1, 
it is possible to bound L 2 [ 0 ,T]^ - norm of the kernel; 

( 2 . 22 ) 


||^i||L2[o,r]2 — 


kl{s, u)dsdu = 

pT 


kl{s, u)dsdu + 


k\{s, u)dsdu 


'0 Jo 
< C 


10 JO 




Jo J u 


'0 Jo 

pT 


< Chi,H2 


pT 




dsdu) 


0 J u 


U 






u 


1 - 2 H 1 


(T - u) 


AH2—^H\ —1 


du 


< < 00 . 

It means that the integral operator KTf{u) = ki{s,u)f{u)du is compact lin¬ 
ear self-adjoint operator from L2[0,T] into L2[0,T] and Fredholm alternative can 
be applied to equation fl2.2ip . To avoid the question concerning eigenvalues and 
eigenfunctions, we produce the following trick. 

It is very easy to see that for any a > 0 


K{ta,sa) = K{t,s)a2~^^^ , dtKHi,H 2 {tci-, sa) = dtK{t, s)a 2 ^^^ 

k{ta,sa) = a^^^~'^^^k{t,s), 


(2.23) 

whence 


ki{ta, sa) = ki{t, s)a'^^‘^ k 

Therefore we can put in equation fl2.2ip s = s'T, u = u'T and hT^z) = hT(Tz), and 
equation fl2.2ip will be reduced to the equivalent form (we omit superscripts) 

/•i 


{uT) 2 ='hrp{u) + / 'hT{s)ki{s,u)ds = 'hT{u) 


(2.24) 


-l-A / hT{s)ki{s,u)ds, 


with A = 2 " 2 ir 2 - 2 Hi ^-2 operator Ki is compact linear self-adjoint operator 

from T2[0, T] into T2[0, T], as it was mentioned above, it has no more than countable 
number of eigenvalues any of them are real numbers, and with only one possible 
condensation point 0. Taking the sequence T„ —)■ cxd in such a way that 

\ _ n^'iH\—2H2^2 

^n- Jn IHi 

will be not an eigenvalue, we get that equation fl2.24p with as upper bound of 
integration has unique solution whence the proof follows. □ 


Now we establish the form of maximum likelihood estimate. 


Theorem 2.4. Let H 2 — Hi > Then the likelihood function has a form 

(2.25) Lt{X, 9) = exp{96H,N{T) - ^9^Sl^{N){T)}, 

and maximum likelihood estimate has a form 


(2.26) 


9{T) 


N{T) 

^HAN){Ty 
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where N(t) = Eo{Xi{t))\^^) is a square integrable Gaussian -martingale, N{T) = 
Jq hT{t)dX{t) with hT{t)t^~^^ G L2[0,T], hT{t) be a unique solution to fl2.9p and 


Proof. We start with fl2.2l) . Consider Gaussian process N{t) = EQ{Xi{t)\'^f). 
Since Xi{t) is ^^rmartingale and C the process iV is a -martingale 
with respect to probability measure Pq- Furthermore, we can present V{t) as 
V{t) = Eo{X‘f{t))\^f) — N‘^{t). Note that Xf{t) — 2 - 2 H i?r mart ingale. Therefore, 
(2.27) 


Po N\t) - 


t 


2-2H 


2-2H 


-G(t)) ) =Po(i?o(X^(t)|l?f)- 




t 


2-2H 


= Po Xl{t) - 


f 


2-2H 


2-2H 


2-2H 

..2-2H 


d: 


X 


df = Eo{Xtis)\^f)- 


2-2H 


= AC 


..2-2H 


2-2H 


V{s] 


and this means that the quadratic variation of the martingale N equals {N){t) = 
j-‘ 2 .- 2 H _ yfif and the likelihood ratio is reduced to 


(2.28) Lt{X, 6) = exp {96h,N{T) - ^9^Sjj^{N){T)y 

so, we get (12.2511 and (I2.26p . Now, taking (12. 9 p into account, we get that 


(iV)(T) = Po(iV2(T)) = Po(y^ hT{t)dX{t)') =Eo[j^ hT{t)d{X,{t) + X2{t)y 

2 r '^ 2 r '^ 

= Po(/ hT{t)dX,{t)'^ +Po(/ hT{t)dX2{t)'^ =^ll hl{t)t^-^^^dt 


rT rT 


hT{u)duKH^^H2{'^i'b)du / hT{s)dsKHi,H2{^i'b)dsdt 


'0 Jt 


I'sAu 


= 7^1 / h‘^{t)t^ -|- / hT{u)hT{s) I dsKHi,H2{s,'t)duEHi,H2{u,t)dtdsdu 


= 7^1 / hT{t)G ‘^^^dt, 


whence the proof follows. 


□ 


In what follows, saying “T cxo” we have in mind that the corresponding property 
holds for any sequence —)■ cxo that has only hnite common points with the sequence 
of eigenvalues of operator Ki. Proof of the following result repeats the proof of the 
corresponding statements from |1] so is omitted. 


Theorem 2.5. The estimator 9 t is unbiased and the corresponding estimation error 
is normal 


9t- 9^ NlO, 


/„ hT{s)s^~'^^^ds 


Now we establish the asymptotic behavior of the estimator. 
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Theorem 2.6. Let H 2 — Hi > j. Estimator 6 t is strongly consistent and 
lim -9? = -. 

/g hQ{u)u2 

Proof. At first we rewrite equation fl2.2ip in the equivalent form, changing u = 
u'T, s = s'T and omitting superscripts: 

(2.29) + f hT{sT)ki{s,u)ds, 

Jo 

or 

(2.30) [ hT{sT)T^^-hi{s,u)ds, 

Jo 

Denote /i = = fihT{uT)T^^~^. Then, taking into account 

fl2.23p . equation fl2.29p can be rewritten as 


(2.31) 
Note that 


1 1 

= -7ifiV(“) + / hf,{s)ki{s,u)ds. 




'0 


{N){T)= j hT{s)s^-‘^^Hs= I hT{s)s2-^^ds = T^-^^^ 

Jo Jo 

Define the operator K 


h^{u)u^^ ^^du 


{Kf){u)= / f{s)ki{s,u)ds 
Jo 

and the scalar product {f,g) = f{s)g{s)ds, f,g & -^ 2 ( 0 , !]• Then equation fl2.3ip 
can be rewritten as 

„2 1 


(2.32) 

Note that 

(2.33) 


+ Kh^{u). 


^•1 rl 


ff 


{KfJ)= {Kf){t)f{t)dt= (/ f{s)ki{t,s)ds)f{t)dt 
Jo Jo Jo 

rsAt 

f{t)t^^~^^‘^f{s)s^^~^^‘^ / dsKHi,H2is,v)dtKH^^H2{t,v)dvdsdt 


n 




= dv dsKHi,H2{s,v)f{s)s^^ / dtKHi,H2{tH)f{t)t^^ ^^^dt>0. 

J 0 J V J V 

Introduce corresponding the first type auxiliary integral equation 
(2.34) = {Kh){u). 

It follows from Lemma [XT] that fl2.34p has the unique solution, say, hg, obviously, not 
depending on /i. The function 5^ = h^ — hQ satisfies two equations + = 0 

and K5fj, + ~ Multiplying the 2nd equation by and integrating, we 

get 
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and it follows from fl2.35j) and fl2.33j) that 5/^11^ < | (ho, 5^)1 < || ho|| ||5/^11, which 
implies that ||(5^|| < ||ho||. Multiplying the 1st equation by ho and integrating we get 

{K6^, ho) + ^o) = 0. 

jJ. 

Note that inequality ||(5^|| < ||ho|| implies that 

Kh;„ho)| < ho)| + ||ho||^ < 2||ho||^ < oo, 


and hence 

(2.36) = I (.5,., Mo) I = l(A'3„fc„)| = i7l,.l(/v.'"o)l ^ 0 

/i 

as T —)■ cx). It means that limr-s-oo Jq h^{u)u^~^^du = ho(u)u^~^^du. Therefore 


T 


2 - 2 H 2 


Ee{9T - 9f = 


-)■ 


hfj_(u)u2 ^^du Jq ho{u)u^ ^^du 


whence the proof follows. 


□ 


Remark 2.7. In outline, our method of proof follows the method of the corresponding 
result from [3j, however Lemma [3.11 is specihc to our case. 


3. Appendix 


We recall some notions from fractional calculus. For the details see na. Fractional 
integrals are dehned as 

(/“+/)(x) = [ f{t){x - t)°-^dt 

r(cr) Ja 

and 

{Ib-f){x) = [ f{t){t - x)°-^dt, 

r(a) Jx 

while fractional derivatives are dehned as 

and 

Fractional differentiation and integration are inverse operators on the appropriate 
functional classes. Also, we shall use the following integration by parts formula for 
fractional derivatives, 

[ {V2+f){x)g{x)dx = f f{x){V^_g){x)dx. 

J a J a 

Lemma 3.1. For any constant C > 0 integral equation 
(3.1) = C{Kh){u), u e (0,1] 

of the 1st kind has the unique solution. 
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Proof. We can present eqnation fl3.1j) in eqnivalent form 


u 


1/2 Hi _ (j h{s)ki{s,u)ds), mG ( 0 , 1 ], 


or 


fsAu 


u ^=C h{s) dsKHi,H 2 {s,v)duKHi,H 2 {u,v)dvds, u e {0,1], 

Jo Jo 

where h{s) = h{s)s^^‘^~^^, or, at last, 

(3.2) h{s)dsKHi,H 2 {s,v)d^duKHi,H 2 {u,v)dv. 

Now, taking into acconnt the transition from eqnation fl2.8p to fl2.9p with the help of 
differentiation, we can perform the inverse operation and get from fl3.2p the following 
eqnivalent eqnation 

(3.3) 

=C{2-2Hi) Kh,,h,(u,v)(^J h{s)dsKH,,Hc{s,v)dsjdv,u€ [0,1], 

The right-hand side of equation fl3.3p can be rewritten as 

PU 

C{2-2Hi) KHi,H 2 {u,v)q{v)dv, 

Jo 

where q{v) = h{s)dsKH.,^^H.^{s,v)ds. At first, solve the equation 

PU 

^2-2Ri^^ / KHi,H 2 {u,v)q{v)dv, 

Jo 

with Cl = C{2 — 2Hi). Taking into account fl2.4p . the latter equation can be 
rewritten in equivalent form 


. 2 - 


U 


^H^=Cif3H2 {u-zY/^-^^z^^-^^{z-vf^-^/^dzq{v)dv, 


or 


y2 ‘iHi ^ I ^H2 


H 2 -H 11 


1I/2-H1 / „1/2-H2 


yi/2-H2^^ - vf^-^/‘^q{v)dvdz. 


'0 


or, at last, 

where C 2 = CiI 3 h 2 ^{‘J>/‘J ~ H 2 ) and 

(3.4) p{z)=z^^-^^f v^l‘^-^^{z-v)^^-^l’^q{v)dv. 

Jo 

It means that 


(3.5) 


p(p) = C,-'(P“^"-(-"-“'))(«) 

= (CjE)//, - l/2))“‘( f (u- t)“'-^'^r-'‘"'dt'j' = 


where C 3 = ^ 1 / 2,3 2 Ri) ^ Furthermore, comparing fl3.4l) and fl3.5lh we get 

that 

^3^1/2-112 ^ f yi/2-R2(^ _ y f 2-3/2^ 
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whence 

= CsiTiH^ - l/2)r(3/2 - f/ 2 ))"'( J (v- 

where C 4 = Cs{2 — 2H2){T{H2 — l/2)r(3/2 — H 2 ))~^. Obviously, q{v) = 
and we arrive at the equation 

(3.6) = [ h{s)dsKH^^H2{s^v)ds. 

J V 

Note that 

dsKH,Ms,v) = /9H.r(3/2 - _ ^)H,-3/2 j j 

so, with the help of integration by parts formula, equation fl3.6p can be rewritten as 

O = j' (■ - 

= iV{H2 - l/2))-^{lfr^^‘^iV^r^^‘^h)-^^-^^){v), 

where ~ The latter equation means that 

= C(i{l - 

h{v) = v^^-^/^h{v) = 

and this solution of equation fl3.ll) is unique. □ 

Acknowledgement 3.2. My thanks to M. Kleptsyna for the statement of the prob¬ 
lem and to M. Kleptsyna and G. Shevchenko for helpful discussions. 
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